Analytical results from the quantum theory of a single-emitter nanolaser. II. 



O 
(N 

a 

1—5 

m 
(N 



C/3 

o 

• I— I 

> . 

Oh. 

q 
o 

• I— I 

Oh- 



> 

in 

o 

m 



X 
J3 



Nikolay V. Larionov^B and Mikhail I. Kolobov^ 
^Department of Theoretical Physics, St-Petersburg State Polytechnic University, 195251, St.-Petersburg, Russia and 

^ Laboratoire de Physique des Lasers, Atonies et Molecules, 
Universite Lille 1, 59655 Villeneuve d'Ascq Cedex, France. 
(Dated: January 25, 2013) 

Using two equivalent approaches, Heisenberg-Langevin and density operator, we investigate the 
properties of nanolaser: an incoherently pumped single two-level system interacting with a single- 
cavity mode of finite finesse. We show that in the case of good-cavity regime the Heisenberg- 
Langevin approach provides the analytical results for linewidth, amplitude fluctuation spectrum, 
and intracavity Mandel Q-parameter. In the bad-cavity regime we estimate the frequency at which 
the peak of relaxation oscillations appears. With the help of the master equation for density operator 
written in terms of coherent states we obtain approximated expression for Glauber-Sudarshan P 
function. This solution can be used in the case of good-cavity regime and allows us to investigate 
in more detail the thresholdless behaviour of the nanolaser. All two approaches are in a very good 
agreement with each other and numerical simulations. 



1. INTRODUCTION 

One of the first theoretical work devoted to possibil- 
ity of realization of a single-emitter laser was published 
by Mu and Savage in Three- and four- level pumped 
emitters placed into lossy cavity and interacted with cav- 
ity mode have been considered there. Some interesting 
effects which do not appear in conventional lasers have 
been revealed, such as squeezing, self-quenching in in- 
coherently pumped laser, deviations from the Schawlow- 
Townes formula for the linewidth. Subsequent publica- 
tions related to single-emitter problem [2|-|9| have proven 
these results and also discovered new effects connected 
with the presence of one emitter: vacuum-Rabi doublet 
in the spectrum 0, 0|) lasing without inversion 0, en- 
tanglement between emitter and field subsystems [7|. In 
most of these works results have been obtained with help 
of different type of numerical methods applied to the mas- 
ter equation for density operator or to the Heisenberg- 
Langevin equations. 

Nowadays a single-emitter laser is realized in exper- 
iments where the role of emitter can play the single 
atom or ion 11 1 or quantum dot [l3|- In this con- 
nection, the purpose of our paper is to provide analytical 
results which can be considered as a tools for experi- 
menter. Thus in our previous Brief Report (l3| with help 
of the Heisenberg-Langevin approach we obtained ana- 
lytical expressions for linewidth, amplitude fluctuation 
spectrum and Mandel Q parameter, which describe the 
behaviour of a single-emitter laser in the case of good 
cavity regime. We also improved the condition required 
for thresholdless regime. In this paper in Sec. 2 we will 
consider and analyze in more detail the linearization pro- 
cedure which we used earlier and extend our calculations 
to bad-cavity regime. In Sec. 3 we will write master equa- 
tion for our nanolaser in the coherent state representation 
and in stationary regime derive approximated expression 
for Glauber-Sudarshan P function. With the help of the 



latter we will prove our previous results, analyze some 
limited cases and threshold behaviour of nanolaser. 



2. HEISENBERG-LANGEVIN APPROACH 
2.1 ModeL Semiclassical theory 

The simplest model of a single-emitter nanolaser is the 
single two-level system placed inside the single-mode cav- 
ity and incoherently pumped to its upper level. Just four 
constants characterize this nanolaser: F - the incoher- 
ent pumping rate, pumping process associated with the 
transition |1) |2); 7/2 - the decay rate of polariza- 
tion due to spontaneous emission to modes other than 
the laser mode; k/2 - the field decay rate in the cavity; 
g - the coupling constant between the field and the two- 
level system. The Heisenberg-Langevin equations for our 
nanolaser are 
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-d{t)+ga{t) + fa{t), 
^{T + j)a{t)+gD{t)d{t) + U{t), 

r-7-(r + 7)^(t) 

2g\a\t)d{t) + a^t)ait)\ + Mt). (1) 



Here a{t) and a''{t) are the photon annihilation and cre- 
ation operators in the cavity mode, a = |1)(2| is the 
operator of polarization of the two-level system, and 
D = |2)(2| — |1)(1| is the operator of inversion between 
the upper level |2) and the lower level |1) of the sys- 
tem. We use the fact that for a single two-level system 
|2)(2| + |1)(1| = I - where I is a unity operator. 

The fa (t), fa [t) and fo {t) are the Langevin noise 
operators which arise through the interaction with the 
heat baths. Using the standard quantum-optical meth- 
ods (see, for example, Refs. [l3|)j we obtain the following 
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nonzero correlation functions of these operators: 

(fait) flit')) ^ K5{t~t'), 

{fl{t)Ut')) = mt-t'). 

{miKt')) - i5{t-t'), (2) 

Cmhit')) = 2^{a{tmt-t'), 

{fD{t)Ut')) = -2T{a{t))5{t-t'), 
{fl{t)fD{t')) = -2T{cj\t))5{t-t'), 
UD{t)fl{t')) = 2^{a^{t))5{t-t'), 

{fD{t)fD{t')) = 2[{r + j)-{r-j){bm5{t^t'). 

We are also interested in the field properties transmit- 
ted outside the cavity through the outcoupling mirror. 
For that we should consider the annihilation operator 
o-out{t) of photons outside the cavity [l5| . The input- 
output transformation is 



aout(i) = v^a(t) - fa{t)/\fn. 



(3) 



In the next Sec. 2.2 we will applied linearization pro- 
cedure to Eq. ([T]) which allows us to investigate the small 
fluctuations of laser field near the dominant stationary 
semiclassical mean value. To find the latter for our 
nanolaser we need to collect a large number of photons 
in the cavity and also provide the coherent interaction of 
each photon with the single emitter. 

In the assumption that above two conditions are sat- 
isfied let us obtain the stationary semiclassical solu- 
tion from Eq. ([T]). For that we need to drop all time 
derivatives dX/dt and average the remained stationary 
equations under assumption of following decorrelations 

^i)a^ = ^-D^ (a), (o'^a) = (o"^) (a)- Suppose that in 
steady state the optical phase is randomly distributed 
between and 27r we take the average value of opera- 
tors with the fixed arbitrary mean value of the phase to 
be equal to zero ipo = 0. Such selection of the phase 
value implies the reality of semiclassical field amplitude 



ao and polarization ctq 



^9 



qq. Thus, the analytical ex- 



pression for semiclassical stationary intracavity intensity 
^0 = I flop is therefore 



/n = 



(r-1) 



(r + l)2 



(4) 



where we introduce new dimensionless parameters: the 
dimensionless pumping rate r = r/7; the dimensionless 
saturation intensity 1^ = "f /k; the dimensionless coupling 
strength c = Ag'^/nj. 

Eq. (01 coincides with one firstly obtained by Mu and 
Savage /q (r) is a parabolic function of pump rate 
which has physical interpretation when c > 8 and when 
the value of the pump rate lying in the domain between 
two points, so called threshold rth and self-quenching 
points, which are given by the following expressions 



where = c/2 — 1 is the point where the stationary 
solution has a maximum I„i = Is (c/S — 1). 

The stationary intensity Iq (r) increases as a pump 
is increased and it starts from the threshold point rth- 
When c ^ 8 then rth ~ 1 + 4/c, what indicates on fi- 
nite value of the threshold or, in other words, there is no 
thresholdless regime in this semiclassical model. The self- 
quenching point rq corresponds to pumping rate where 
the atomic polarization is rapidly damped to zero due to 
emitter trapping in the excited state, what leads to the 
damping of the field. 

From the expression for the maximum of intensity /,„ it 
follows that large mean number of photons in the cavity 
can be achieved when clg ^ 8 (we assumed that c ^ 8). 
The latter condition can be written in other form g ^ k. 
Thus the semiclassical regime takes place when the pho- 
tons lifetime in the cavity is so long that each photon 
is provided by coherent interaction with the single emit- 
ter (above mentioned condition). The condition clg ^ 8 
also gives us possibility to collect a large mean number of 
photons in two opposite cases, namely in good - /s 3> 1 
and bad-cavity regimes Ig <^ 1. 



2.2 Linearization around the stationary solution 

To linearize the Heisenberg-Langevin equations ([T]) we 
assume that all operators can be presented as a sum of 
dominant classical term Xq and a "small" operator val- 



ued fluctuation SX (see, for example, Refs. |l6l. Il7j). 

X ^Xo + 5X. (6) 

After linearization with respect to 5X we obtain the fol- 
lowing equations for the fluctuations 



-5d{t) = - -5a{t) + g5a{t) + /Jt), 



dt 

d „ , , 

d 



dt 



6D{t) 



- ^{T + -f)6a{t) + g\ao6D{t) + Do5d{t) 

+ Lit): 

~ {r + j)5D{t)-2g\<Jo{Sa{t)+5a^{t)) 

+ ao{6a{t)+Sa^{t)\ + fnit). (7) 



Before solving Eqs. ([7|) we first split operators into Her- 
mitian "real" and "imaginary" parts as 



6a{t) = 6u{t) +iSi){t), 
Sa{t) = Sfiit) +iSri{t), 
fit) = t{t)+iA{t). 



(8) 



rth = r„ 



-yjl - 8/c, rq = rra + 



(5) In this way the linearized equations of motion separate 
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into two independent blocks 
d 



dt 
d 

di 



6u{t) 



5u{t)+gSij.{t) + Y.a(t), 
(r + -i)5pL{t) +g\aQ5D{t) + Da5u{t) 



2 



j^5D{t) = - (r + -f)5D{t) - 4g \aa5u{t) + a^5ji{t) 
+ Id it); 



(9) 



"imaginary" parts is a result of our phase selection in 
the derivation of semiclassical solution (see Eq. (|4]) and 
text above). 

The block for three real parts 6u, Sfi, 5D are related 
to the intensity fluctuations via 51 — 2ao6u. The two 
imaginary parts Sv, Sfj can be associated with phase fluc- 
tuations. 

Eqs. ([HI [T0| can be resolved by means of Fourier trans- 
formation. For that we need to perform the fluctuations 
as 



dt 
d 

dt 



Suit) 



Siy(t) + gSflit) + Aait), 

{T + -f)Sm+9DoSHt) + Kit). 



(10) 



We want to note that the independence of "real" and 



5X{t)= 5X{n)eicp{~int) 



dn 

2^' 



ill) 



which gives us the linear algebraic equations which can 
be resolved by means of Cramer's rule. The result for the 
Fourier-transformed amplitude 5u and phase 50 quadra- 
ture components are 



in 



{in-{T + j + K))/2){in - (F + 7)) + V|ao|2 



gA^{n)-Aain)iin-(T + -f)/2) 



in 



in - (F + 7 + k)/2 



4Kg2|ao| 



r 



(12) 
(13) 



where 

A{n) 
c{n) 



~m - (F + j)/2){in - (F + 7)) + Vlaop], 
g{in-{T + ^)), B = -g^ao. (14) 



Eqs. and from Eq. pTjl the Fourier-transformed 
quadrature components Su, Si> are (5-function correlated. 
Using the input-output transformation Eq. ^ we find 
the amplitude and phase fluctuation spectra 



As it follows from the correlation function definition 



{6uoutin)duont{.^i')) = -sin + n') 
{Si>out{n)Si)ont{n')) = ^s{n + n') 



4 172(02 -d)2+ 72 (el]2_/)2 
■2 ^2^cir + 1) + {r + 1)^] 



4 n'-{n^ + j^r + i + i/is)yA) 



(15) 
(16) 



r 



where 



-f^[c{r - 1) - {r + if + 2cr], 



b = -f^ir~l)[3{r + lf + 2c 
d = 



2 c/,(r-l) + (r + l) 
^ 2L 



c(r + l)(r-5)], 
3(r + 1) + 



(17) 



Now we will use these results to analyze two different 
regimes of our nanolaser, namely good- and bad-cavity 
regimes. 
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2.3 Good-cavity regime 

In the case of good-cavity regime Is ^ 1 the ob- 
tained expression for the amphtude fluctuation spectrum 
Eq. becomes 



1 



S{r,c) 



1 + [n/noy 



(18) 

The first term in the square brackets corresponds to 
standard quantum hmit (SQL), the second term is a 
Lorentzian with width f2o given by 



1 - 



(r + l)2 



c{r - 1) 



(19) 



and this term describes the additionally exes above the 
SQL or reduction of the noise below SQL depending of 
the sign of the strength 5(r, c) 



S{r,c) = 



{r - 1) [3(r -t- 1)3 -f- 2c2 - c(r + l)(r - 5)] 
[c(r-l)-(r+l)2]2 



(20) 

We have investigated this expression in order to find 
out a possibility of the noise reduction below the SQL in 
the laser light outside the cavity. Unfortunately, this kind 
of nonclassical phenomenon are very limited. We have 
found that in the region of pump parameter ?',„ < r < 
the strength S'(r, c) is always positive and increasing func- 
tion, resulting in the excess noise in Eq. (|18p . The small 
negative values of S{r, c) are observed when the pump 
parameter is centered around r = c/5 (c > 200), i.e. 
lies in the region rth r < Vm- Probably, these negative 
values are the result of the well-known antibunching phe- 
nomenon for a single-emitter. However, this antibunch- 
ing effect is strongly attenuated due to the effect of ac- 
cumulation of many photons inside the cavity during the 
long (and random) photon lifetime. 

The Mandel Q-parameter can be expressed through 
the amplitude quadrature component 5u as Q = 4, ■ 
J {S u (fl) 5 u{~ri))dfl/ 271 — 1 and the result is 



1 



1 - 



(r + l)2 



c(r - 1) 



(21) 



We see that the strength S{r,c) is also characterized 
this quantity. On the Fig. la there are some graphics 
for Mandel Q-parameter for different coupling constant 
c when Jg ^ 1 . For the pump rate lying between thresh- 
old and quenching points our theoretical results are in 
a very good agreement with the exact solution of the 
master equation for the density matrix (see next Sec. 3). 
The main divergency appears near the threshold and 
quenching points where the linear theory does not work. 
Two fluctuations peaks of Mandel Q-parameter are well 
known [Hj : first of them associated with laser turn on and 
second broad peak associated with laser turn off. 



Another quantity of interest is the low-frequency 
asymptotic version of the fluctuation spectrum of the 
phase quadrature component. Using the Eq. (fT6|) we ob- 
tain 



1 -|-c/(r + 1) 



(22) 



The low-frequency divergence in this spectrum as l/fl^ 
is manifesting for the phase diffusion process [l^l- This 
spectrum also defines the linewidth Aiy of our nanolaser 
through {{Si'out f}n/Io 



2 \ r + 1 



(23) 



where Az^st ^ i^/ (2/o) is the Schawlow-Townes linewidth 
of a conventional incoherently pumped laser. The ob- 
tained formula for the linewidth shows the deviation from 
the Schawlow-Townes result and this coincides with the 
numerical simulations performed by Clemens et al 

At the end of this section we want to discuss the le- 
gality of above linearization procedure. Let us explain it 
in terms of adiabatic elimination of the emitter variables 
in Eq. ([1]). Indeed, in the case when Is ^ 1 one can ob- 
tain the Heisenberg-Langevin equation only for the field 
operator a{t). The semiclassical solution of the latter 
coincides with Eq. The following linearization with 
respect to fluctuation dd{t) gives us equation which con- 
tains the products like da{t)fa{t) {a ^ a, D). Because 
of presence just one emitter they can not be neglected in 
frame of linearization procedure. However, these prod- 
ucts appear in the equation as terms multiplied by small 
factor 2K/(r — 7) < I/Iq ^ 1, what makes the order of 
the products similar to second order in fluctuation Sa{t). 
Finally, the fully linearized equation for a(t) gives the 
same results as was obtained in this section. 



2.4 Bad-cavity regime 

The considered single-emitter laser relates to so called 
high (3 lasers - the lasers with high fraction /? of spon- 
taneous emission into the lasing mode (see expression 
for /3 in Sec. 3, Eq. (j44|l ). Most of present work on 
the intensity noise in high /3 lasers has been applied to 
standard semiconductor laser diodes, embedding quan- 
tum wells as gain material. These models predict that 
if /3 is small, the peak-to-valley difference of the relax- 
ation oscillations is very large and it decreases with in- 
creasing /3 1^, 2^- However, all these models are based 
on standard rate equations used in quantum well lasers 
and very few studies have been carried out in this fleld 
when switching from quantum well lasers to quantum dot 
lasers. Some recent research was focused on the turn-on 
dynamics and relaxation oscillation in standard quantum 
dots lasers |2l| but has not been extended to unconven- 
tional high /3 nanolasers. 
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Figure 1: a) The Mandel Q-parameter (|2ip vs the pump pa- 
rameter r. Dashed Une shows results obtained with help of 
function Po JMI- 1) c = 20; 2) c = 30; 3) c = 40. For all 
curves Is = 20. b) The spectrum of amplitude fluctuations. 
I) r = rm = 199, 2) r = 3. For all curves Is = 0.2, c = 400. 



Here we use results obtained in Sec. 2.2 to analyze be- 
haviour of our nanolaser in the case of bad-cavity regime 
Is 1- As was mentioned above, for that we only need 
to preserve the following condition cIs ^ 8. 

If the pump rate r is close to the maximum point rm 
the amplitude fluctuations spectrum is described by the 
same Eqs. (|18[[^^ as in a good-cavity regime (see Fig. lb, 
curve 1). The main difference is the magnitude of the 
width Qq ()19p for the amplitude fluctuations spectrum. 

Other situation takes place when the pump rate r lies 
not far from the threshold point rth ■ In this case the high 
peak in the amplitude fluctuations spectrum appears (see 
Fig. lb, curve 2). The physical origin of this peak is the 
relaxation oscillation phenomenon. Because of lineariza- 
tion procedure break down in the vicinity of threshold 
rth , we can not correctly define width and height of this 
peak, but we can estimate the frequency of the relaxation 
oscillations flosc/n = Is\/c{r — l)/2. 



3. MASTER EQUATION APPROACH 

This section are devoted to investigation of the 
nanolaser properties with help of master equation for 
density matrix. This approach have been used by many 



authors. Thus in [l[ the master equation have been 
reduced to system of first order, ordinary differential 
equations and analyzed numerically. In [5|, |9| more ef- 
ficient quantum trajectory algorithm have been devel- 
oped for numerical evaluating of master equation. The 
works 01 are devoted to analysis of the master equa- 
tion written in terms of coherent states for field. In 
light of our interest, we want to briefly discuss one of 
them, namely Authors have worked with the sys- 
tem of equations for Glauber-Sudarshan P function and 
additional quasi-probabilities. Besides numerical simula- 
tions they have also found approximated expression for 
P function which can be used in the case of good-cavity 
regime. Here we follow the nearest approach as in Q, 
but we manage to derive isolated stationary equation for 
P function. A detailed analysis of the latter allows us to 
obtain approximated solution, which has wider region of 
application than that in Q . 



3.1 Coherent state representation of master equation 

The master equation for our nanolaser is (see, for ex- 
ample, 0) 



dp 
at 



V,P 



(2 apa^ 



a' ap — pa' a 



(2a'' pa — aa^ p — paa'') , 



(24) 



where interaction between the cavity mode and the sin- 
gle two-level system is given by the Jaynes-Cummings 
Hamiltonian 



V = ihg (a^CT — (T^a) . 



(25) 



Finally, we want to obtain equation for P function, for 
that we need to rewrite Eq. ([M)) in terms of coherent 
states for field \z), \z*) and in projections of the density 
operator on the two- level system states \1) and |2). Using 
well known rules (see, for example, Ref. [22|) 



ap 
pa 

t 



pa 



a'p -> 



zp [z, z 
_ _d 
dz 

z p{z,z ) 

* „ A 

dz 



p{z,z*) , 



/5(z,z*) , 



(26) 



where z, z* are the complex variables and intro- 
duce the following quasi-probabilities pik (z, z*) = 
{i\piz,z*)\k) = p^k («,fc = 1,2), D = P22{Z,Z*) - 
pii{z,z*) and P function P ~ pii{z,z*) + p22{z-,z*), 
we obtain from Eq. the system of partial differential 
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equations 

dP d (n 



dt dz\2 

(r-7)P-(r + 7)^ + ^(f^^ + 



d /K 
d / K 



{t^zP - gp2i ] + TjTT ( ■7;z*P - 5P12 



dt 



9P21 



9p; 



'21 



(r + 7) 



dt 



9 



(9 / K 

(9z 



(27) 



The first equation in this system is written in special 
form which can be associated with conservation law for 
the quasi-probability dP/dt = divJ, where J is known 
as probability current density. The additional quasi- 
probabilities D and p2i have a clear physical meaning. 
Thus, the mean values for inversion and polarization of 
the two-level system are expressed through D and p2i as 
D)=fD (z, z*) (Pz, (a) = / P21 [z, z*) (Pz. 



3.2 Stationary solution 

In stationary regime the probability current density J 
can be equated with zero (see, for example, [23, 24|) what 



gives the following relations p2i = (P12)* = {i^/'^9)zP. 
The latter allows us to eliminate p2i, P12 from the sta- 
tionary form of the system (|27p and obtain two coupled 
equations for P and D. In the new variables / = |z| and 
if {z = \/l exp (*</j)) these equations read 



/P - I [(r -l)P-{r + l)D] = ~I{P + D), (28) 



^p {r + 1 + 1/ Is) 
c 

2 81 cIs y ^ I 



ID+\{P + D) 



where we introduced above mentioned dimensionless pa- 
rameters Is, c, r. 

There is no phasing influence on our nanolaser, e.g. 
- external field with fixed phase or coherent pumping. 
Therefore expected stationary solution does not depend 
on phase and we can write P and D as only function 
of / = |z|^. This is also following from the second 
Eq. (|29p : the term proportional to imaginary unit should 
be equated to zero because of reality of functions P and 
D. 

Let us derive isolated equation for the P function in 
the stationary regime. Using Eqs. (^511^5)1 we can express 
function D in terms of P 



D 



{r + -21 /h\cl^ dl 



d 



-P 



2 [3 - /, (r + 1 + c)] 
c/? 



r 



IP+{r-l + \/Is)P 



(30) 



Now put this relation into the first or the second equation 
in the system Eqs. ([28l [29 |) . In such a way we get isolated 
equation for P function 

P" {I)+p{I)P' {I)+q{I)P{I)^0, 
p{I) (aio + flu/ + ai2p) I {aQ2l^ + aosP) , 
q{I) = (a20 + 021/ + a22/^) / (ao2/^ + aoa/^) , (31) 

where the prime indicates derivative with respect to vari- 
able / = \zf and functions p{I), q{I) depend on physical 
parameters Is, c, r via aik (see Appendix). 

The obtained equation is the second order differential 
equation with polynomial coefficients. To use this equa- 
tion we need to define boundary P{0), P(oo) or "ini- 
tial" P{0), P'{0) conditions. If we equate the variable 
/ to zero in Eq. ((3T|l we obtain the following relation 
P'(0) = — a2o/aioP (0). This relation together with nor- 
malization of the quasi-probability P gives a comfortable 
way to define the "initial" condition for numerical simu- 
lations. 



In the case of good-cavity regime functions p{I) and 
q{I) in Eq. ([3T|) contain large parameters like cIs, 1^ or 
Ig (see Appendix) . Therefore, we can try to find approxi- 
mated solution of Eq. (I3T|) using the perturbation method 
(see, for example, Ref. [25|)- The point of this method 
is representation of the unknown function P as iterative 
series P = X)i^o where A is a large parameter. 
A substitution of the latter series into Eq. (|5T|) with sub- 
sequent equating to zero of terms with same powers 1/A 
gives possibility to find the different orders of approxi- 
mations P„. 

To isolate the correct single large parameter let us 
rewrite p(I) and q{I) in the following form 

p{I) = ai2 (/ - I-i) {I - I+i) //' (/ - loo) , (32) 
qil) = 022 (/ - I-2) (I - 1+2) 11^ [I - loo) , (33) 

where roots are 



iil2aa±\l (0^1/012)^ - 4a,o/ai2/2, 



'00 



-ao2- 
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In the conditions Is ^ 1, c > 8 the magnitudes of the 
roots 7-1-2, loo approximately equal each other and 
the roots /_2 have the similar order. So, it is natural 
to save all I±i and extract the large parameter from the 
constants ai2, 022 in Eqs. p21 I33p . This parameter is 
therefore A = cIs (= a22 ~ —012) and Eq. pip can be 
written as follows 

A-ip" (/) + p{I)P' (/) + q{I)P (I) - 0, (34) 

where p{I) = p(/)/a22, qil) = q{I)/a22- Here we want 
to note that we do not neglect the term r/c in relation 
o-vil 0-12 = (7 — 3/s(r + 1) — 2c/s) /2cls^ because of pump 



rate r can be comparable or bigger than c (for example 
near to maximum Vm or to quenching Vq points). 

Thus, in zero-order approximation (first term in the 
series P = X]^o "^" ) '^^ neglect the term in Eq. ([34]) 
proportional to 1/A what implies first-order differential 
equation for function Pq 

ai2 (/-/-i)(/-/+i)P^(/) 
+ 022 (/ - I-2) {I - 1+2) Po (I) = 0. (35) 

Po is then found from Eq. ([55)) : 



iVo (1 - I/I-l)^' (1 - I/I+lY' CXp (-f2|/) , if / < /_! 



Po(/) = <^ ' ' -12V' ^^^- .^ (36) 

\ 0, otherwise 

_ 022 (/-I - /-2) (^-1 - ^+2) „ a22 {I+i - I-2) (I+i - 1+2) 

Jl — 77 r ^ ' /2 — 



ai2 (/_i-/+i) ' ai2 (/-i-Z-i 



where Nq is normalizing constant. The roots and con- 
sequently the powers /i, /2 in Eq. (|36p have complex 
structure and will be estimated below in some special 
cases. 

To understand why our solution is restricted we note 
that the roots I±i of function p{I) in Eq. ([3T|) are known 
as turning points, where oscillating behaviour of P (/) 
can change on exponential. A careful examination of 
Eq. ([5T|) shows that in the conditions I, ^ 1, c > 8 on 
the /-domain ranging from to turning point /_i the 
P (!) is positive and nonoscillating function. Moreover, 
its essential changing occurs just on latter domain (it 
is not for every set of parameters r, Is, c, see below). 
Otherwise, when / > /_i then P (/) becomes oscillating 
function. The obtained solution Eq. ([36)) can not describe 
any oscillations, moreover it possesses imaginary values 
when / > therefore we have restricted definition 

domain for our solution by the value /_i and call them 
boundary root. 

To reveal physical meaning of the boundary root we 
have simplified expression for /_i and found maximum 
of function Pq (/) on the domain (0,/_i). When Is ^ 1, 
c > 8 we have managed to estimate the latter quan- 
tities with help of results obtained in Sec. 2.3: /_i w 
Ist/ {2 + 3r/c) > Io[l + Q {r,c)] and the point where 
Po (/) reaches his maximum is / = /_2 ~ Iq- Thus, 
P function has a maximum at the point corresponding to 
the semiclassical intensity Iq Eq. ([4]) and the boundary 
root is at a distance IoQ{r,c) from it, where Q{r,c) is 
the Mandel Q-parameter Eq. ([^ . 

In the Fig. 2 we plot the P function calculated using 
Eq. ()36p (dashed line) and using numerical simulations of 



the equation pip (solid gray line). The different curves 
correspond to different values of pump rate r while c 
and Is are fixed. We see a very good agreement with 
numerical results. 

The main divergency appears for value of pump rate 
discovered in Sec. 2.3, namely for r = c/5, where the 
Mandel Q-parameter Q {r, c) ([2T)) has a minimum (see 
Eq. ([20)) and text below) . In the right wing of P function 
(Fig. 2b) the small oscillations appear when c becomes 
comparable with Is, but still our solution is good. With 
growing of c the oscillations become more pronounced 
and the maximum of P function goes out of domain 
(0, /_i) and our solution p6p becomes inapplicable in 
vicinity of r = c/5. It is interesting to note that critical 
situation occurs when c ^ 200. For latter values of c and 
for arbitrary saturation intensity Is the maximum of P 
function goes out of the domain (0, /_i). As mentioned 
above, the P function acquires oscillating character for 
/ > what probably indicates on the nonclassical ef- 
fect associated with photon antibunching. 

Now let us analyze the obtained solution ([36)) in some 
special cases. 



Far below threshold and far above quenching 

At first we consider situation when the pump rate lies 
far below threshold r <^ rth- In this case the average 
value of / = |z|^ is small and the following inequalities 
« // (/,r/2) < 1, J//+1 « I /{Is (r + 1) /2) « 1 
are satisfied. Thus, the factor in front of exponential in 
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Figure 2: a) P function vs I — \z\^ for different pump param- 
eter r. 1) r = 8, 2) r = 16, 3) r = 24, 4) r = 36, 5) r = 44. 
For all curves Is = 100, c = 50; b) Oscillation behaviour of 
the P function. Is = c = 100, r — c/5 = 20. Dashed line - 
solution (|36|) . solid gray line - numerical simulation, dotted 
line - solution from Ref. 0. 



Eq. l[36|) can be approximated as follows 



{i-i/i^,f{i~i/u,y^ 



= exp [/i In (1 - 
sa exp[ 

The Pq is therefore 
^o(/) = 



///_i)]cxp [/2ln(l-///+i)] 
-(/i//_i + /2//+i)/]. (37) 



1 



aio/O'20, 



■ exp 



(aio/a2o) 



(38) 



The obtained expression for P function corresponds to 
thermal distribution with the following mean number of 
photons in the cavity (n) = aio/a2o ^ f- 

For the values of pump rate lying far above quench- 
ing point r ^ Tq we can make the same manipulation 
with Eq. ([36| as above. Using inequalities I/I-i ~ 
//(I«r/ (2 -I- 3r/c)) < 1, I/I+i « //(/sr/2) < 1 which 
occur in this case, we get the same thermal distribution 
as Eq. (j38|) , but the mean number of photons possesses 
another behaviour {fi) = 010/020 ^ 1/r. 

In the two limiting situations, when r — > or r — )■ 
00, the corresponding mean number of photons converge 
to zero and the distribution Eq. gives Dirac delta- 



function Pq (J) = (5 (/), what coincides with the thermal 
distribution with zero temperature. 



Near to maximum point Vm 

Here we consider situation when the pump rate 
is in the vicinity of r„j, where the semiclassical so- 
lution Iq Eq. (|4]) reaches his maximum. In this 
case the powers in Eq. p6p are simplified as /i w 

(3 (r + 1)^ + 2c^-c{r + 1) (r - 5)) / (3 (r + 1) + 2c)^ 
/2 ~ and the root /_2 becomes similar to Iq. Simple 
algebra leads to following expression for Pq 



Po (A/) 



No 



{I- AI/Afcxp{AI) 
0, 



if A/ A 
otherwise 



(39) 



where we introduced new variable A/ = I — Iq and con- 
stants A = fi/a, a = 2c/ (3(r + l) + 2c). For chosen 
values of pump rate the average value of A/ is small and 
the inequality A// A ^ 1 is satisfied. Thus for the factor 
in front of exponential in Eq. p9p we can write 



(1 - A//A)' 



exp [Aln(l - A//A)] 
exp [-A/ - A/V2A] 



what implies Gaussian P function 

(I 



Po (/) = No exp 



2A/a 



(40) 



(41) 



Here we removed the restriction on Po- Indeed, the dis- 
tribution (|4ip never really sees the boundary and / can 
be taken to run from to 00. 

Using the obtained Eq. (j4l]) we can estimate the av- 
erage value of photons in the cavity (h) and its variance 
{{Ahf) 



/ IPo{T)dI^Io, 
Jq 



(42) 



{{Any) = / iI-{h)rPo{I)dI 

Jo 

« /o[l + A/(a/o)]. (43) 

A substitution of the explicit form of A and a into the 
expression for variance gives ((Afi)^) = Io[l + Q {r,c)], 
where Q (r, c) is the Mandel Q-parameter obtained from 
linear theory (|2ip . 

In the weak coupling regime c < Is the obtained 
Eq. (|4T|l works good for values of pump rate lying in 
whole semiclassical region rth < r < Vg. In the strong 
coupling regime c > Ig the symmetrical Gaussian distri- 
bution (j4ip is not sufficient and it becomes more pro- 
nounced near to mentioned value of pump rate r = c/5. 
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Threshold 



In our previous report |13| we considered nanolaser 
behaviour around the semiclassical threshold rth (0) 
and specified the condition required for thresholdless 
regime. Here we continue our research using the obtained 
Eq. pel). 

The behaviour of P function in the threshold point r^h 
for three different regimes 1) c ^ /s, 2) c w /s, 3) c ^ /s 
is shown in Fig. 3a. To realize these regimes we fix the 
large coupling constant c = 100 and set the different 
values of the saturation intensity: Ig = 600, 60, 6 (we 
choose such values of /, for better resolution of P func- 
tion behaviour). In the weak coupling regime (c <C /«) 
the P function has a typical plateau (it is marked by 
solid straight line, curve 1), which indicates transition 
to lasing: from thermal to Gaussian type distribution. 
As saturation intensity is decreased and the strong cou- 
pling regime [c^ Is) occurs the maximum of P function 
moves from zero value of variable / and the semiclassical 
threshold behaviour disappears (curve 3). 

Above we obtained the maximum point for Pq (|36p 
and it was the root /_2. To define the threshold value 
of pump rate we need to solve equation /_2 (r) = 0, i.e. 
find such value of r when the maximum of the P func- 
tion is in the point / = /_2 (r) = 0. The approximate 
solution isf=l-|-4/c — 2/Ig. When c <C then the 
last term can be neglected and f equal to semiclassical 
threshold nh ~ I + i/c > 1 ([5]). When c > /s then 
f becomes smaller than unity what indicates the disap- 
pearance of the semiclassical threshold and transition to 
thresholdless regime (see the vanishing of narrow peak 
in the behaviour of Mandel Q parameter in our previews 
report [l3|, Fig. 2). The dynamics of f in above three 
considered regimes is f sa 1.04, 1, 0.66 (remind that 
rth = 1.04). 

At the end of this section we want to correct the mis- 
print occurred in our previous report (see formula (13) 
in [13 )■ The valid expression for fraction /3 of sponta- 
neous emission into the laser mode is 



/3 = c/ [(c+l) + /,('^ + l)] 



(44) 



We should note that this misprint did not affect on all 
our results and discussions. 



Comparing with other authors 

In introduction to this section we mentioned that close 
approach based on quasi-probabilities have been consid- 
ered by two authors Karlovich and Kilin in They 
have worked with first order differential equation system 
for function P and additional quasi-probabilities. By ne- 
glecting of small parameters in this system they have ob- 
tained integrable equations and as a result the analytical 



P{/) 

0,4 
0,3 
0,2 
0,1 
0,0 



P(') R 

"■3 ■ \ 5 



a) 



10 



20 



30 



/ 40 



0,2 



0,1 



\ 



/ 



\ 

X 



1 b) 

3 // \ \ 

\ 1 1 



1 



V / / \ 



/ X. 



2 4 



10 12 14 / 



Figure 3: P function vs J = |2| . a) Transition to the thresh- 
oldless regime. 1) = 600, 2) Is = 60, 3) 7s = 6. For all 
curves c = 100, r = rth ~ 1.04. b) Parameters was taken 
as in Ref. 0: 1) c = 2000, = 5, r = 6, 2) c = 1428.57, 
Is ^ 7, r = 4.28, 3) c = 1000, L = 10, r = 3, 4) c = 500, 
Is = 20, r = 1.5, 5) c = 285.7, Is = 35, r = 0.86. Dashed line 
- solution (|36[l. solid gray line - numerical simulation, dotted 
line - solution from Ref. 



formula for P function. The structure of their solution 
has an identical form with our Eq. (|36p (see formula (9) 
in Q), but this solution works good in the region near to 
semiclassical threshold rth or quenching points r^. For 
the pump rate lying around the maximum point rm it 
gives uncorrect result for width of function P and as a 
consequence the uncorrect variance of the photon distri- 
bution. 

In the Fig. 3b we plot some graphics of quasi- 
probability P for parameters taken from [6| . It is clear to 
see that our Eg. ( [5^ gives the same result as a solution 
(9) obtained in [6| and all of them are in a good agree- 
ment with numerical simulations. The Fig. 2a, b demon- 
strate that for the pump rate lying around the maximum 
point r^n our Eq. (j36p (dashed line) has an advantage 
with (9) Q (dotted line). The latter indicates well the 
maximum of true P function, but the width is uncorrect. 

The derivation of a single equation for P function (PT|) 
and accurate extraction of large parameters gives us pos- 
sibility to obtained the solution which has more wide re- 
gion of application. 
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4. SUMMARY 
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In this paper we have studied physical properties of a 
single-emitter laser. The problem have been investigated 
in terms of the Heisenberg-Langevin equations and in 
terms of the master equation for density matrix. In both 
approaches we have provided analytical results which are 
summarized below. 

In the case of good-cavity regime with help of the 
Heisenberg-Langevin approach we have obtained analyt- 
ical expressions for linewidth Eq. p3|. amplitude fluc- 
tuation spectrum Eq. (|18|) and Mandel Q parameter 
Eq. ((2T|) . These results work good in the semiclassical 
region of pumping rate rth < r < rq. According to 
nanolaser behaviour the latter region can be split into two 
subregions. In the first subregion r < Tm the nanolaser 
behaviour is similar to that of conventional laser. Also 
in this subregion for r centered around r = c/5 and for 
large coupling constant (c ^ 200) we have discovered the 
small noise reduction below the SQL in the laser light 
outside the cavity and small negative values in Q (r, c) 
(j2ip . The master equation approach confirms this non- 
classical phenomenon: P function manifests oscillating 
behaviour. Probably, these negative values are the result 
of the well-known antibunching phenomenon for a single- 
emitter. In the second subregion r > the excess noise 
is always takes place. The P function in this subregion 
does not have any nonclassical features. 

In the bad-cavity regime we have observed two differ- 
ent situations: when the pump rate r is close to its maxi- 
mum value Tm then laser generates as in the good-cavity 
regime; when r is close to threshold rth then high peak in 
the amplitude fluctuations spectrum appears which indi- 
cates on the relaxation oscillation phenomenon. 

With the help of the master equation for density matrix 
written in terms of coherent states we have managed to 
derive the stationary equation for the Glauber-Sudarshan 
P function (|3H) . A detailed analysis of this equation al- 
lowed us to obtain approximated solution Eq. p6p , which 
works good when clg ^ 1. We have analyzed Eq. ([55]) 
in some special cases. In the good-cavity regime in the 
semiclassical region rth < r < Tq Eq. ()36[1 can be written 
as Gaussian distribution with mean value coincides with 
the semiclassical intracavity intensity /q (|3]) and with 
width IqQ (?', c) (HU. For the values of pump rate lying 
far below threshold and far above self- quenching points 
Eq. f36|) corresponds to thermal distribution. 

Using Eq. ([55]) we have also obtained approximated 
expression for threshold pump rate f = 1 + 4/c — 2//^. In 
the weak-coupling regime {la c) the last term can be 
neglected what implies the semiclassical threshold f = 
rth ~ 1 + 4/c. When strong-coupling regime {Is <C c) 
occurs then f < rth what indicates on the transition to 
thresholdless regime. 
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APPENDIX 



The constants aik in Eq. ([21 



ao2 
ao3 

aio 
On 

0-12 
0-20 



a2i 

0-22 



1 _ (r + 1) 

2^2' 
1, 

ocr _ 3 cr(r + l) 
^ 4^ ^ 4 

9 _ 6(r-H) + c ^^ i{r+lY + c{Ar + 2) 
4^2+^ 4 
7 _^ 3(r + 1) + 2c 

2 " 2 ' 

3 ll(r + l) + 3c 2 'i{r + lf + 2c 
2 " 4 ^ " 2 



,c(r + l) 



(r-1)- 



{r + lf 



= c/. 



(45) 
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